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I. INTRODUCTION 



Quantum degenerate exciton-polaritons in a semi- 
conductor microcavity have recently attracted growing 
interest. Such a system consists of quantum wells sand- 
wiched between two distributed Bragg reflectors, where 
excitons in the quantum wells are coupled with pho- 
tons confined in the reflectors. The coupling of the ex- 
citons with photons makes the exciton-polariton mass 
much smaller than the electron mass, which enables 
the system to maintain long-range coherence. In recent 
years, experimental studies on this system have shown 
rapid progress^: for example, Bogoliubov excitations^ 
and quantized vortices^^ have been observed, and vari- 
ous quantum fluid dynamics have been demonstrated^^ — 
In the present paper, we focus on a system in which po- 
laritons are continuously pumped by near-resonant laser 
beams. ^^^^^ The loss of polaritons with a short lifetime is 
balanced by the continuous pumping, thereby realizing a 
nonequilibrium open system. 

Homogeneous steady states far from equilibrium can 
spontaneously break the translational symmetry, and 
they can exhibit patterns or structures, called dissipative 
structures. This kind of structure formation occurs in 
a variety of systems, such as Rayleigh-Benard convec- 
tion cell^ in a fluid heated from below, and Belousov- 
Zhabotinsky patterns and Turing patterns in chem- 
ical reaction-diffusion systems. However, few studies 
have been performed on dissipative structures in quan- 
tum mechanical systems. ^^^^^ The polariton superfluid in 
a microcavity is a suitable system for realizing quantum 
dissipative structures, since polaritons can be injected 
into the system continuously and coherently with rela- 
tive ease.^ The continuous injection and dissipation of 
polaritons is expected to drive the formation of quantum 
dissipative structures. 

In the present paper, we show that various steady pat- 
terns of quantized vortices are formed in an exciton- 
polariton superfluid in a semiconductor microcavity 
pumped by near- resonant continuous- wave laser beams. 
We consider a situation in which two spatial regions are 
pumped and a polariton superfluid flows into the inter- 
space between the two pumped regions, as illustrated in 
Fig. [U (a). When the interval d between the two regions 
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FIG. 1: (color online) (a) Schematic illustration of the sys- 
tem. The quasi-two-dimensional system lies in the x-y plane 
and is pumped by laser beams applied from the z direc- 
tion. This coherently excites the polaritons with zero mo- 
mentum. The pumped regions are restricted to \y\ > d/2 
and the pumped polariton superfluid flows into the interspace 
1^1 < d/2. (b) Dispersion relation of the lower and upper po- 
laritons. The pump energy and momentum correspond to the 
circle. 



is small, the steady state has continuous translational 
symmetry (with respect to x in Fig. [D (a)). When d ex- 
ceeds a critical value, an instability sets in and the con- 
tinuous translational symmetry is spontaneously broken. 
The system develops into another steady state, in which 
quantized vortices are aligned periodically. These pat- 
terns of quantized vortices are spontaneously organized 
in nonequilibrium steady states and are peculiar to the 
open system with loss and pump. We also show that the 
pattern depends on d and the phase difference between 
the two pumped regions. In finite systems, vortices es- 
cape from the edges of the pumped regions, which can 
be plugged by an appropriate potential. 

This paper is organized as follows. Section Ull provides 
the formulation of the problem. Section IIIII numerically 
demonstrates the formation of various patterns. Sec- 
tion |IV] presents the conclusions to this study. 



II. FORMULATION OF THE PROBLEM 

We used the mean-field theory to study the dynamics 
of an exciton-polariton superfluid. The energy of a pho- 
ton in a quasi-two-dimensional microcavity can be ap- 
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proximated as TiuJq +p^/(2Mc), where p is the in-plane 
momentum and Mq is the effective mass of a photon. The 
effective mass of an exciton is much larger than Mc, and 
we approximated the energy of an exciton as a constant, 
Tiujy^. The excitons in the quantum wells and the pho- 
tons in the microcavity couple with the Rabi frequency 
r^R. The macroscopic wave functions for the photons ?/^c 
and the excitons ?/^x are assumed to obey the coupled 
nonlinear Schrodinger equations given by^— 



ifi 



dt 



dt 



2Mc 



(la) 



where V{r) is the external potential for the photons, g 
is an exciton-exciton interaction coefficient, and 7c (x) 
is the decay constant for the photons (excitons). The 
last term on the right-hand side of Eq. ([Taj) describes the 
coherent pumping of photons by external laser beams, 
where /Cp, cjp, and F{r) are the in-plane wave number, 
frequency, and profile of the pumping, respectively. 

Diagonalizing the first and second terms on the right- 
hand side of Eq. ([1]) for the mode with wave number 
we obtain the energies of the upper and lower polaritons 
as 



where 



(2) 

hk'^ /{2Mc)' The dispersion relation in 
Eq. (|2]) is depicted in Fig. [TJb). The difference between 
the pump and lower polariton frequencies is defined as 



S = Up — a;_(0). 



(3) 



For simplicity, in the following calculations, we restrict 
ourselves to the case of uJq = ujy^ = 0, 7c = 7x = 7, 
and kp = 0. The parameters are fixed as follows: Mq = 
2 X 10~^me with rrie being the electron mass, fil^R = 5 
meV, g = 0.01 meV/im^, 7"^ = 10 ps, and S = 0.76 
meV. 



III. NUMERICAL RESULTS 

A. Steady states with continuous translational 
symmetry 

We first consider an infinite system without an external 
potential (^ = 0). We use a pump profile as illustrated 
in Fig. [U^a), given by 



Fo {y < -d/2), 
F{r) ={ {\y\< d/2), 
Foe'^ {y>d/2), 
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FIG. 2: (color online) Steady states having continuous trans- 
lational symmetry with respect to x, which are obtained by 
solving Eq. ([T]) with pumping as in Eq. (|4]). The phase differ- 
ence in the pump is = in (a) and (f) — 71 in (b). The main 
panels show the minimum value of iV^c {y)f of the steady state 
(corresponding to the density dips as shown in the inset). For 
d < dmax, there exist stable (red solid) and unstable (green 
dashed) steady states. The density and velocity distributions 
of the photonic states corresponding to the circles are shown 
in the insets, where the red solid and green dashed curves 
indicate the stable and unstable steady states. The vertical 
lines in the insets indicate zLd/2. 



where Fq = 5.7 meV, (j) is the phase difference between 
the two pump beams, and d is the width of the interspace 
between the two pumped regions. For this value of Fq, 
the pumped regions are excited into the stable uppermost 
branch of the bist ability curve. 

In the limit of d ^ with (/) = 0, the whole system 
is uniformly pumped into the stable state '0c,x(^,^) oc 
g-zwp^ When d is increased, the density at \y\ ^ d/2 de- 
creases and the continuous translational symmetry with 
respect to x is maintained if d is not very large. We 
sought such steady states that is uniform in x by solv- 
ing Eq. ^ using the Newton-Raphson method. We also 
checked the stability of the steady states by using Bogoli- 
ubov analysis. If there is a Bogoliubov frequency that has 
a positive imaginary part, the state is dynamically un- 
stable. Figure [2] shows the d-dependence of the steady 
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(a) d = 24 |im 
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(b) d = 30 |im 





states. For ^ = 0, the steady states have even-parity 
density profiles with two minima, where the velocity is 
largest, as shown in the inset in Fig. [2fa) (the exciton 
wave function ?/^x is always qualitatively the same as the 
photon wave function ipc^ and thus in the following dis- 
cussion, we will show only i/jc)- There are two steady 
states for each d; one is stable (red solid curve) and 
the other is dynamically unstable (green dashed curve). 
The dynamically unstable modes grow exponentially and 
break the continuous translational symmetry with re- 
spect to X. The density dips in the unstable state are 
deeper than those of the stable state. These minimum 
values of the density are shown in the main panels of 
Fig. [21 As d is increased, the stable and unstable branches 
meet at d = (imax and the steady states disappear for 
d > <imax, which is a manifestation of the saddle- node 
bifurcation. Figure Efb) shows the case of ^ = tt. The d- 
dependence is similar to that in Fig. [2fa) except that the 
density and velocity profiles are asymmetric, as shown in 
the inset of Fig. [2fb). There are two degenerate states, 
since Eq. ([1]) is symmetric with respect to y ^ —y and 
i^c,x — '0c,x for (/) = TT (only one state is plotted in 
the inset). 

From Fig. O we expect the following scenario: When 
the initial value of d is small and then d is increased 
adiabatically, the system follows the stable steady states 
(red solid curve in the main panels of Fig. [2]). When d ex- 
ceeds (imax, the steady states disappear, which causes the 
breaking of the continuous translational symmetry with 
respect to x. The system then goes to another nonequi- 
librium steady state with broken symmetry, if it exists. 
In the next subsection, we will show that these states 
have a variety of patterns with quantized vortices. This 
scenario of pattern formation is similar to that in the 
dissipative structures'^. 

The disappearance of the steady states for a large d 
can be understood as follows. The steady states have 
the factor of e~'^^p^ and therefore have the "energy" of 
hujp. Since the kinetic energy of the superflow must be 
smaller than this energy ficjp, the flow velocity has an 
upper bound fmax- The loss of polaritons per unit time at 
y < \d\/2 is ^ jpd^ where p is the polariton density. The 
loss must be replenished by injection ~ pv, i.e., jpd ~ pv. 
Thus, d has an upper bound dmax ^ ^max/7- We have 
numerically confirmed that (imax is roughly proportional 
to 7-^ 



1 




(c) d = 33 |im 








(d) d = 40 urn 
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FIG. 3: (color online) (a)-(d) Density (iV^d^, left panels) 
and phase (arg-^c, right panels) profiles of the steady states 
for (a) d = 24 /xm, (b) d = 30 /xm, (c) d = 33 //m, and a 
snapshot of an unsteady state for (d) d = 40 /xm. The unit 
of the density is 350 /im~^. The flow fields in the dashed 
squares in (b) and (c) are shown in (e) and (f), respectively. 
The two pumped regions have the same phase, = 0, which 
corresponds to Fig.[2l^a). See Supplemental Material^^ for the 
dynamics, in which d is changed from 24 fim to 40 fim. 



B. Symmetry breaking and pattern formations 

Figures [3](a)-[3]^c) show the stable steady states for 
= 0. The steady state shown in Fig. [3l^a) is uniform in 
X, which corresponds to the stable branch for d < (imax 
in Fig. [2fa). As d increases and exceeds (imax, the stable 
branch in Fig. [2fa) disappears. As a consequence, the 
continuous translational symmetry with respect to x is 
spontaneously broken and quantized vortex pairs are cre- 
ated; these align periodically, as shown in Fig.[3fb). The 



vortex state shown in Fig. [H^b) is a stable steady state, 
in which vortex dipoles (clockwise and counterclockwise 
vortex pairs) are the building blocks of the pattern: two 
vortex dipoles form a pair that has parity symmetry with 
respect to y. Although a vortex dipole propagates at a 
constant velocity in a uniform system, the vortex dipoles 
in Fig. [3fb) remain at rest. This is because the vortex 
dipoles are prevented from propagating by the counter- 
flow from the pumped regions, as shown in Fig.[3je). The 
polaritons pumped in the region oi y > d/2 {y < —d/2) 
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FIG. 4: (color online) (a) Density (|'0c|^) and phase (arg'^c) 
profiles of the steady state for d = 24 /xm, where the phase 
difference between the two pumped regions is = tt. The 
unit of density is 350 /im~^. (b) The flow field in the dashed 
square in (a). See Supplemental Material^^ for the dynamics 
in which is changed from to tt. 



flow in the direction of —y {-\-y) and sustain the vortex 
dipoles located in y > {y < 0). These vortex dipoles 
would propagate in the -\-y {—y) direction if they were in 
a uniform system without superflow. 

As d is increased further, the steady vortex structure 
in Fig. [3l^b) makes a transition to another vortex struc- 
ture, as shown in Fig. [3fc). The vortex pairs aligned 
along y are not vortex dipoles but twin vortices in 
which two vortices have the same circulation. Such twin 
vortices with clockwise and counterclockwise circulations 
alternate along y — 0, resulting in a structure similar to 
the vortex street reported in Ref. [23|. Although two vor- 
tices with the same circulation rotate around one another 
in a uniform system, the twin vortices in Fig. ^c) are at 
rest. This can be understood from the Magnus force on 
the vortices, which is due to the density distribution^^ 
in Fig. [3fc): a vortex with circulation n experiences an 
effective force proportional to —k x Vp, where p is the 
polariton density. For example, in a uniform system, the 
twin vortices magnified in the right panel of Fig. [3ff) 
would rotate around one another in the counterclockwise 
direction. The density gradient in the directions of the 
red (light gray) arrows exert forces in the directions of 
the blue (dark gray) arrows, which holds the rotation of 
the twin vortices. 

The continuous translational symmetry with respect to 
X, as shown in Fig. [3l^a), thus changes to discrete trans- 
lational symmetry due to the vortex pattern formation, 
as shown in Figs, ^h) and[3fc). For a larger (i, there 
is no steady state, and the vortex pattern dynamically 
changes in a random manner, which has no symmetry (a 
snapshot of the dynamics is shown in Fig. [3fd)). 

Figure m shows the steady state for ^ = tt, which cor- 




FIG. 5: (color online) Density profiles (|'0c|^) for the finite 
pumped regions given by Eq. (|5]). (a)-(d) Snapshots of the 
dynamics for (f) — tv and d — 24 fim. The origin for the 
time is chosen appropriately. The arrows trace the vortex 
motion. The dashed circles indicate a created vortex pair, 
(e), (f) Steady states with the plug potentials given by Eq. (|6]) 
with Vo — 380 meV for (e) = tt and d = 30 /xm, and (f) 
= and d = 24 /im. The plug potentials are marked by 
the dashed squares. The unit of the density is 350 /im~^. See 
Supplemental Material^^ for the dynamics in (a)-(d). 



responds to the case of Fig. [2fb). For the width d = 24 
/im, there is no steady state that is uniform in x, and a 
vortex pattern emerges as shown in Fig. HJ The clock- 
wise and counterclockwise vortices are aligned alternately 
along ^ ^ 0, whose flow pattern (Fig. HJb)) is similar to 
that in Fig. [3ff). For a classical incompressible inviscid 
fluid, such an arrangement of point vortices is unstable. 
(The only stable arrangement of this kind of vortex rows 
is the one found by von Karman.^^) The vortex pattern 
in Fig. HJis stabilized by the superflow from the pumped 
regions \y\ > d/2 toward y = 0. The density distribu- 
tion around the vortices may also be responsible for the 
stabilization. 

The vortex patterns in Figs. [3fb), [3l^c), and S] are pe- 
riodic in X, and their periods are :^ 28 /im, 31 /im, and 
34 /im, respectively. The periods of the vortex patterns 
are determined by the geometry of the system, and they 
are of the order of the interval d between the pumped 
regions. 



C. Finite systems 

Thus far, we have considered only infinite systems, in 
which the steady states have continuous or discrete trans- 
lational symmetry with respect to x. In this section, we 
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examine pumped regions of finite size, as follows: 

{Fo {\x\ < 60/im, -40/im <y< -d/2), 
Foe'^f" (|x| < 60/im, d/2 < y < 40/im), 
(otherwise), 

(5) 

where Fq = 5.7 meV is the same as that in Figs. [3] and HI 
Figures [5fa)-[5fd) show the dynamics for the parameters 
(p = 7T and d = 24 /im, where the origin of time is chosen 
appropriately. For these values of d and 0, the infinite 
pumped system has the steady state with a vortex pat- 
tern, as shown in Fig. |4j In the finite system, however, 
the polaritons in the space between the pumped regions 
flow out from the edges at \x\ = 60 /im. As a conse- 
quence, vortices also flow out, as indicated by the arrows 
in Figs. [5ja)-[5fd). As the vortices flow out, they sep- 
arate from each other, and new vortex- ant ivortex pairs 
are created between them, as shown by the dashed circles 
in Figs.[5fc) and[5fd). Such separations and creations of 
vortices are repeated, and the system does not reach a 
steady state. 

To stabilize the system and obtain steady states, we 
add a plug potential as 

y(^) ^fVo (50/im < \x\ < 60/im, \y\ < -d/2), 
^ ^ \0 (otherwise), ^ ^ 

which prevents the polaritons from flowing out. Fig- 
ures [5fe) and[5ff) show the steady states for (/) = tt and 
^ = 0, where the regions of V = Vo are enclosed by the 
dashed squares. The clockwise and counterclockwise vor- 
tices alternately align in Fig. [5je) and the vortex dipoles 
make pairs in Fig. [5ff), which are similar to Figs. |4]and 
[3jb), respectively. Thus, the steady states can also be 
obtained for finite systems by using a plug potential that 
confines polaritons to a finite region. The steady-state 
patterns are suitable for time-integrated imaging in ex- 
periments. 



ductor microcavity; this is a nonequilibrium open sys- 
tem. We considered the situation in which polaritons 
with k = are pumped into two regions and flow into 
the space between them, as illustrated in Fig. [TJ When 
the distance d between the two pumped regions is small, 
there exists a stable steady state with continuous transla- 
tional symmetry with respect to X J as shown in Fig.[3fa). 
As d is increased, this steady state disappears (Fig. [2]), 
and the system goes into other steady states with peri- 
odic patterns of quantized vortices (Figs. [3fb), [3fc), and 
SJ. Such patterns are specific to a nonequilibrium sys- 
tem with loss and pump, and they can be recognized as 
dissipative structures. We have shown that steady-state 
patterns can also be observed in a finite system with an 
appropriate potential (Fig. [5]). 

The structure formations shown in the present paper 
are peculiar to nonequilibrium steady states and should 
be distinguished from those in equilibrium systems, such 
as crystals. For example, the triangular lattice structure 
of quantized vortices in a rotating superfluid cannot be 
regarded as a dissipative structure, since it is the ground 
state in the rotating frame of reference. 

Note added. Very recently, a preprint^^ has appeared, 
which reports experimental and numerical observations 
of self-arranged vortex-antivortex pairs in a polariton 
condensate. 
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